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Abstract 

A simple relativistic model of heavy-quark-light-quark mesons is proposed. 
In an expansion in inverse powers of the heavy quark mass we find that all 
zeroth and first order heavy quark symmetry relations are satisfied. The main 
results are: 

• the difference between the meson mass and the heavy quark mass plays a 
significant role even at zeroth order; 

• the slope of the Isgur-Wise function at the zero recoil point is typically 
less than —1; 

• the first order correction to the pseudoscalar decay constant is large and 
negative; 

• the four universal functions describing the first order corrections to the 
semileptonic decay form factors are small; 

• these latter corrections are quite insensitive to the choice of model param- 
eters, and in particular to the effects of hyperfine mass splitting. 



0. Introduction 

It has recently been realized that major simplifications occur in the theoretical 
treatment of weak decays of mesons containing a heavy quark. ||T| New symmetries 
appear in an expansion in inverse powers of the heavy quark mass. For example 
at lowest order the six form factors for the decays of pseudoscalar to pseudoscalar 
and vector mesons become related to a single universal form factor (the Isgur-Wise 
function). The latter depends only on the dynamics of the light degrees of freedom. 
In addition the masses of pseudoscalar and vector mesons with the same heavy/light 
flavor content are equal at this order, as are their decay constants. 

However in the real world it is important to know the magnitude of the l/mq cor- 
rections. It is expected that the b quark is probably heavy enough for the corrections 
to be small, but it is less clear for the charm quark. It is also becoming apparent that 
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the size of the corrections depends strongly on what quantity is being calculated. We 
will attempt to shed more light on these questions. 

We will present a model based on simple finite quark loop graphs. It is relativistic 
and incorporates the effects of hadronic recoil in a natural way. We demonstrate 
explicitly that all heavy quark symmetry relations among the semileptonic decay 
form factors at zeroth and first order in l/mg are satisfied by the model. The model 
also allows calculations to all orders in I/uiq. 

We consider our approach to be complementary to QCD sum rules. The latter 
approach relies on a certain representation of the hadronic contribution to a 3-point 
function. The basic input for our model is a representation of a Bethe-Salpeter 
amplitude for the heavy-quark-light- quark meson. As such, it may be improved by 
future QCD-based numerical computations of Bethe-Salpeter amplitudes. Our model 
is also complementary to nonrelativistic quark models, 0, ^ and it provides the first 
indication of how the I/tuq corrections are affected by hyperfine mass splitting. (Note 
that hyperfine mass splitting effects are not included in the ISGW model. 0]) 

Our model highlights the importance of the difference between a meson mass and 
the mass of the associated heavy quark. This mass difference is needed to obtain 
consistent results at zeroth order in l/mq {e.g. the Isgur-Wise function). We will 
show that ignoring this difference leads to very different results. And correspondingly 
the l/mq corrections to this mass difference must be properly incorporated into the 
l/rriQ corrections of other quantities. 

Another result is a confirmation that the l/mq correction to the B meson de- 
cay constant fs is very large and negative. This was first suggested by lattice 
calculations, and it has been noted in two dimensional calculations]^ and sum 
rule calculations ^. This is to be contrasted with the corrections to the weak de- 
cay form factors, which we find to be relatively small. The latter observation also 
agrees quahtatively with sum rule results, although the actual numerical values of 
various l/mq corrections differ. 

We test the sensitivity of our results to the parameters of the model, and in 
particular we find that the form factor corrections are quite insensitive. We feel that 
the quantitative differences between our model and sum rules give an indication as 
to the true uncertainty in present theoretical determinations of these quantities. 

1. Definition of the model 

Before beginning we make the standard definitions of physical quantities. The 
meson decay constants are given by 

{0\A^\P{p))=zfpp, and {0\V^\V{p)) = fyMye^. (1) 

Our normalization is such that = q'^^i'^zQ. The various form factors for semilep- 
tonic decay are defined by 



{V2{V2) \V,\ PiK)) = yM^Mv^/ly(c^) W^r^2X, (3) 
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{V2{v,)\A,\Pr{v{)) = 

where the w's are meson velocities, and u = vi ■ V2- 

We choose to represent transition amphtudes by diagrams with heavy mesons 
attached to a loop involving heavy and light quarks. The essential nonperturbative 
physics of QCD is that contributions with large momentum flowing through the light 
quark line are suppressed. We will model this physics by including factors in the 
vertices which damp the loop integrals when the light quark momentum k is larger 
than some scale A set by QCD. This has two desirable effects: the integrals are made 
finite and we may consider k small in comparison to a meson momentum. This latter 
fact is the crucial ingredient that gives rise to the correct heavy quark symmetry 
relations. 

The effective vertices between a light quark, a heavy quark and a pseudoscalar or 
vector meson are taken to be 

^ > 75 and \ 7^, (5) 



where k is the momentum of the light quark. The light and heavy quarks are assigned 
standard propagators with masses and rriQ respectively. With these choices, and 
with n an integer or half-integer, standard methods involving Feynman parameters 
may be used in the one loop computations. 

We may now consider expanding various quantities to first order in A/mg. 

A±=Afl + (±/i-^7)— 1 (6) 



Ak 



) (7) 




/w = a-— 1 +(,„•— I (9) 



The P *-^'V symmetry in d^-j^) at zeroth order is a direct consequence of the model, as 
shown in Appendix A. The same is true of the nontrivial zeroth order mq dependence, 
and in particular the standard scaling fpy oc vtlq'^I'^ . The zeroth order constants A, 
c and a will be completely determined in terms of Kjmq and n without further input. 

The first order constants Bpy, dpy and bpy depend in addition on the values of 
g and h. These latter constants model the effect of hyperfine mass splitting. Through 
their effect on dpy they determine how the pseudoscalar and vector masses approach 
a common value in the heavy quark limit. To first order in l/mq we will find that g 
and h dependence cancels out of certain physical quantities. As with A, c and a, we 
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stress that the quantities -Bpy, dpy and bpy are not parameters of the model and 
will be determined by A/rrig, n, g and h. 

For certain physical quantities we will need to know rUg and A separately. It 
is reasonable that the appropriate effective light quark mass should be of order its 
constituent mass, in common with most successful quark models. As for A we note 
the following relation which follows from (§) 

"^^^i^ I zeroth order ^ 2^' ('^^) 

The first equality is the standard definition of A appearing in the literature. When 
necessary to choose A, it will chosen so as to give reasonable values of A. The 
parameter n determines the form and extent of the damping due to the meson vertex 
factors, and one of our goals will be to study the sensitivity of l/mg corrections to 
the choice of vertex factors. 

The various weak decay form factors may also be expanded, 

hi{(jj) = ai^{uj) + 6chi{(jj)— + 6bhi{uj)— (11) 

rric rrib 

with = ay = = c^As = 1 and a_ = = 0. $,{uj) is the Isgur-Wise function 
and we will calculate it in terms of A/rriq and n. In QCD it has been shown that 
numerous relations exist between the first order corrections. |l^ We demonstrate 
in Appendix B that these relations are satisfied by our model for any values of the 
parameters. This consistency with heavy quark symmetry at first order in A/mq is 
a nontrivial test for any quark model of heavy mesons. 

We will present numerical results for the first order corrections and their first 
derivatives, all at = 1. We do not treat the perturbative QCD corrections; they 
have been calculated elsewhere ||10[| and should be added to our results. 



2. Zeroth order results 

We first consider the results of our model at zeroth order in A/tuq. Details of 
these calculations may be found in Appendix A. Computation of c involves finding 
the zero of the meson two-point function as given in (|33|). Some examples of these 
"mass functions" are displayed in Fig. 1. When we set the location of the zeros equal 
to = ttIq + cArriQ we find the following values of c: 



A/ nig 


Cn=l 


Cn=3/2 


Cn=2 


2 


1.5797 


1.3949 


1.3006 


4 


1.4635 


1.1976 


1.0545 



(12) 



We see that increasing n decreases M. A is constrained to be greater than nig, since 
otherwise we find no sensible zero of the mass function. We will see later that the 
two choices shown in ([l2|) , A/rUg = 2 and 4, imply a reasonable range for A. 

We display in ([T3|) the results for the zeroth order parameters A and a appearing 
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in Zpy and fpy respectively and given by (|36| ) and (pSl). 







^n=3/2 




a'n=3/2 


fl?i=2 


2 


0.9181 


0.9184 


0.9059 


0.25 


0.17 


4 


1.1461 


1.1358 


1.1077 


0.29 


0.19 



(13) 



We do not include in (p!3f ) results for a„=i because in this case the dependence of jp 
on A is not of the form of (^ (it is logarithmic), while fy diverges. For n = 'i/2 
mi, = 4.8 GeV and A = 500 MeV and 1 GeV we find fs = fp* = 40 MeV and 132 
MeV respectively at zeroth order. The strong A dependence refiects fpy oc 

^3/2 

at 

zeroth order, fpy become smaller for larger n. 

The Isgur-Wise function of our model is given by (^2[). It depends only on 
the index n and the ratio a = rriq/A. Let us see the effect of ignoring the zeroth 
order mass difference A. Putting c = in (|^ allows one to perform the relevant 
integrations explicitly, yielding a linear combination 



c=0. 



UJ 



/i(a)ei(a;) + /2(a)6(^) 



/i(«)ei(l) + /2(«)6(l) 
of two functions independent of n and a, 



(14) 



I + UJ 

The coefficients are given by 



2 InUu + \JuS^ - 1 

and ^2(1^) = - 



a \ 2n-l 

a \ 71 



a 



and 



d 



(15) 



2n-l 



t=l 



a ln(t/a^ 



t 



(16) 



4=1 



The expression (^) for ^(co") with nonzero c must be evaluated numerically. We 
plot the result in Fig. 2 together with the c = form in (0). The difference is 
quite significant; this demonstrates the importance of retaining the zeroth order mass 
difference A. A simple quark loop model was recently described by deRafael and 
TaronpT[| in which they obtained ^(cu) = ^2(^^)- But that model, like our model with 
c = 0, does not correctly take into account a nonzero A. 

It is interesting to see how the with nonzero c depends on a. We plot in Fig. 
3 the first and second derivatives of ^{u) at = 1, for the case n = 3/2. We see that 
there is a range of a over which they are quite insensitive to the value of a . Our two 
choices a = 1/2 and a = 1/4 more or less fall within the region of insensitivity. 

Although we cannot obtain ^(o;) analytically we have checked numerically that it 
has the correct behavior at large u, namely it tends to zero. We have also computed 
the derivatives of ^{u) at ci; = 1 : 



A 


,n 


e(l) 


e'(l) 


2 


3! 


r'd) 

4! 


2, 


3 
2 


1 


-1.241 


1.370 


-1.459 


1.536 


4, 


3 
2 


1 


-1.118 


1.047 


-0.917 


0.777 


2, 


1 


1 


-1.359 


1.691 


-2.094 


2.631 



:i7) 
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Our values for ^'(1) are consistent with the present experimental determinations. |]I2 
(Note that these determinations assume some functional form for ^(cj) which is not 
identical to ours.) We see that our model, like most other models, does not satisfy 
the constraints on ,^ (1) and ^ (1) recently argued for in [Tl[]. It is amusing that the 
constraints are satisfied by the inconsistent c = model. 

3. Order- A/mQ corrections 

The purpose of this section is to find the ©(A/mg) corrections and to explore 
their dependence on the parameters, n, A, g, and h. For (A/mq,n) we take the three 
sets of values (2, |), (4, |), and (2, 1). And for each quantity we calculate we will 
indicate explicitly the g and h dependence. 

We remark that just as the zeroth order A dependence of the meson masses {i.e. c) 
plays an important role in the zeroth order results, the next to leading A dependence 
of the meson masses {i.e. dpy) plays an important role in the first order corrections. 
The meson masses enter the loop calculations via the on-shell external momenta 

= Mpy. A dependence from meson masses also originates in the factor My in the 



definition of the vector decay constant (|I]) and in the factors Mp^Mp^ and y Mp-^^My^ 
in the definition of the form factors (Q-j^). 

Our results for the corrected masses Mpy and the corrected vertex normalizations 
Zpy {i.e. the parameters dpy and Bpy) are 



A/rriq, n 




dp 


dv 




2 ^ 

^' 2 

2,1 


0.379 - 
0.385 - 
0.461 - 


{h + g)Q.U<o 
(/i + ^)0.886 
(/i + ^)0.801 


0.335 + (/i - 
0.262 + (/i- 
0.365 + (Zi- 


^)0.646 
^)0.886 
^)0.801 


h./mq, n 




Bp 







:i8) 



' 2 

4 ^ 

^' 2 

2,1 



(19) 



0.266 -{h + g)l.m 0.301 + {h - g)l.m 
0.193 ~ {h + ^)0.989 0.274 + {h - ^)0.989 
0.404 - {h + c/)1.341 0.535 + {h - c/)1.341 

The values of the correction coefficients for the decay constants as defined in 
(I) are displayed in (|20|) . Note in particular the large negative values of hp for the 
pseudoscalar. (From our parameter fit at the end we find that both g and h are 
positive.) The implication is that the first order corrections to fp are nearly of the 
same order (with opposite sign) as the zeroth order values. We conclude that the 
A/mg expansion for fp is breaking down. 



A/niq, n 




hp 


bv 




2 ^ 

^' 2 


-2.92 - 


{h + g)l.S2 


-1.37+ (/i- 


^)1.82 


4 ^ 

^' 2 


-3.82- 


{h + g)l.^2 


-1.43 + (/i- 


g)1.52 



(20) 



We now turn to the first order corrections to the weak decay form factors. It is 
standard PIl to write the first order terms in (ITIT 



as 



'yi{uj)^{uj) = 5chi{uj)— + 5hhi{uj) — . 

rric mb 



(21) 
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In QCD the first order terms may be expressed in terms of tlie universal functions 
Xi,2,3{^) and ^3(0;) in addition to ^(ct;).[|5|, |13| These relations are reproduced in Ap- 
pendix B. Also in Appendix B we demonstrate that these relations are satisfied by our 
model for any values of the parameters. We also verify in detail Luke's Theorem ||13|| 
which reads = X3(l) = 0. 

At = 1 we find the following results in which we temporarily set g = h = 0. 







X2(l) 


X3(l) 


6(1) 


2 - 

^' 2 





-0.153 





-0.002 


^' 2 





-0.132 





-0.008 


2,1 





-0.148 





-0.005 


rriq ' 


X'i(l) 


X'2(l) 




e;(i) 


2 ^ 
^' 2 


-0.368 


0.212 


-0.058 


-0.097 


^' 2 


-0.307 


0.161 


-0.057 


-0.116 


2,1 


-0.419 


0.218 


-0.073 


-0.101 



(22) 



(23) 



When g and h are nonzero the results in ( p2|) do not change. The g and 
dependence of various preceding quantities has canceled out. Only two of the first 
derivatives are changed by the following amounts. 

5xA^)=x{h-2g) (24) 

5x,{l) = -\xh (25) 

X = 0.193,0.003,0.155 for {A/mg,n) = (2, |), (4, |), and (2,1) respectively. In fact 
X2(^) and ^3(0;) are independent of g and h for any u, as demonstrated in Appendix 
B. 

We may translate these results into values for the 7^(1) 's. At this stage we are 
finally forced to make a choice for rriq. We choose rriq = 250 MeV as a representative 
value of an effective constituent quark mass in our loops. Using = 4.8 GeV, 
rric = 1.44 GeV and g = h = we obtain the following results for 7i(l) and 7j(l) 
expressed as percentages of unity. 



A,n 


7+(l) 


741) 


7^(1) 


7A,(1) 


7A.(1) 


7^3(1) 


0.5 GeV, 1 





-8.5 


15.8 





-19.5 


11.1 


1.0 GeV, 1 





-14.8 


27.1 





-31.6 


17.5 


0.5 GeV, 1 





-9.7 


17.9 





-21.8 


12.3 


A,n 


7;(i) 


I'M) 


7^(1) 


7^,(1) 


7^,(1) 


7;,(i) 


0.5 GeV, 1 


-12.9 


-1.7 


-8.4 


-1.2 


8.4 


-2.2 


1.0 GeV, 1 


-20.8 


-3.6 


-11.3 


0.7 


14.0 


0.5 


0.5 GeV, 1 


-20.1 


-2.1 


-11.2 


-3.2 


9.3 


-3.9 



(26) 



(27) 



Again the only effect of nonzero g and h are in the following first derivatives. 

^7^(1) = ^7^,(1) = 57^3(1) = -2Aa; + '-^] (28) 
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6j'4l) = -2Ax{h + g)(— + —) (29) 

X takes the same values as above. 

An important point is that the zero recoil values of the four universal functions 
in (^) and the corresponding first order corrections to the form factors in ( p6l) are 
independent of hyperfine mass splitting, as modeled by g and h. Another is the weak 
dependence of the results in ( p2|) and (^) on the parameters n and A. The results in 
(p6|) and (p7|) show more variation and reflect the fact that they are proportional to 
A. Corresponding to the three rows of these tables we have A 350, 600, 400 MeV 
respectively. 

In our discussion thus far we have been more concerned with the sensitivity of 
the results to the parameters rather than trying to find an optimal set of parameters. 
The latter could be accomplished in the following way. We may calculate a A_|_ for 
each of B* and D* and a A_ for each of B and D by fitting the zeros of the full mass 
functions to the physical meson masses. By fitting these four A-i-'s to the first order 
form in (P) we find, for n = 1, A = 667 MeV, g = 0.057, h = 0.37, and for n = 3/2, 
A = 818 MeV, g = 0.047, h = 0.32. In both cases A = 504 MeV which coincides with 
a sum rule estimate p|, ^. 

Our results for the corrections are somewhat different from a sum rule calculation p| 



which gives ^3(1) = 0.33 and X2(l) = and from an improved sum rule calculation ||14 
which gives X2(l) = —0.038. Note that in QCD X2(l) measures the effect of a 
chromo- magnetic moment operator insertion. Our values for 7i(l) and 7j(l) may 
be compared directly with the corresponding results in 

In this paper we have presented the zeroth and first order model results for various 
quantities of interest in the meson heavy quark effective theory. The same model 
allows the calculation of physical quantities to all orders in the I/ttlq expansion. This 



comparison, to pursued elsewhere, ||15|| should shed further light on the usefulness of 



the heavy quark quark expansion in B and D meson physics. 

Appendix A 

We treat in this Appendix the zeroth order pieces of 

1. Masses and normalizations 

We begin with the pseudoscalar and vector mass functions Tpy^p"^) defined by 
setting the self-energy graphs equal to iTp and —ig^uTy + ■ ■ ■ respectively, where the 
ellipsis denotes p^Pu terms. The light quark momentum k may be chosen to be the 
same as the loop momentum, and the heavy quark momentum is then vMpy + k 
where v"^ = 1. In the heavy quark limit we may ignore k compared with Mpy and 
use to show that the heavy quark propagator becomes (after scaling k Ak) 

i J^' . (30) 

A —2k ■ V — cpy 

We have temporarily allowed the constants cpy describing the lowest order difference 
between quark and meson mass to be different for pseudoscalar and vector. Our first 
task is to show that they are in fact equal to a common constant c. 
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The relevant traces are 





(^ + 1) 1 (^ + «) (31) 



where a = iriq/A. The translation is k —>■ k — xv where a; is a Feynman parameter; 
the integral linear in k then vanishes. We may anticommute the resulting v leftwards 
past the gamma matrices since terms proportional to do not contribute in the 
vector case. The results are 



a; + «)Tr ^ (^ + 1) 1 \=d , (32) 




where ^ = —k ■ f + a is equivalent with x + a. We thus find 
2 f d^k 4zn,Z|ryA2-4« ^ 

Tpy{Mpy) - J (2^)4 [_p + l]2n [_p + «2] [_2fc . ^; - Cpy] ' ^^^^ 

We are regarding these mass functions as the full two-point functions of the 
mesons; i.e. there are no tree level contributions. Thus the meson masses are de- 
fined by the locations of the zeros of the real parts of these mass functions: 

Re Tpy{Mj,y) = 0. (34) 

Notice that imaginary parts are present in any quantity which is evaluated on meson 
mass shell. This is because the difference A between the meson mass and the heavy 
quark mass is greater than the light quark mass rriq, and thus the meson is above the 
threshold for two free quarks. But note that our A is consistent with that of other 
approaches. P, R| We will not consider these imaginary parts further. We also find 



and then drop an overall, physically irrelevant minus sign. [|T5 

We regard cp and cy as variables whose value is fixed by ( |33D and (^4|). It is then 
obvious that they are equal to a common value c because the functions Tpy{Mpy) 
differ at most by multiplicative factors. This demonstrates that Mp = My at zeroth 
order (see (|)). 

The mass functions are also required to satisfy the normalization condition 

VpyiM^) = 1. (35) 

This immediately implies that the normalization factors Zpy are equal to one another 
in the heavy quark limit. The value of the zeroth order constant A defined in (0) is 
thus fixed by 

^ ^ dI^Ld^ ' ^^^^ 

where Dfc = —k'^ + 1, Da = — fc^ + a^, D = —2k-v — c. Here and in the following we will 
adopt the convention that an overall factor AiUcA^"' and / d^k/ {27:)^ are understood 
whenever the denominator is written explicitly as a product of D's. 
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By the Ward identity the above determination of A is equivalent to the normal- 
ization of the Isgur-Wise function. This will be verified explicitly below. 

2. Decay Constants 

We again choose the light quark momentum to coincide with the loop momentum. 
The traces relevant to the decay constants as defined by (0) at lowest order are 

^ { T ) "* " M X } <^ ^ - { T ) ^ ^' { X } 

= fl I (37) 



Use of the lowest order pieces of the expansions immediately shows that the 

decay constants are equal at lowest order and are given by 

^3/2 _y^-2n^ 

fp,v = X (38) 

3. Form factors 

We turn to the zeroth order results for the meson form factors defined by (@-@). 
The relevant traces are 

Trj J^^ |(^2 + l)r^(^i + l)75(^ + a) (39) 



u 



The translation is k ^ k — xvi — yv2 with = 1, where x and y are Feynman 
parameters. The term linear in k vanishes, and Vi and f 2 may be anticommuted to 
the left and right respectively to yield the standard traces 

{x + y + a)Tr | ] (^2 + l)r^(?^i + 1)75 (40) 



V 



Use of the lowest order pieces of the expansions (|g-@) and evaluation of the standard 
traces immediately yields the relations 

hj^{uj) = hy{uj) = hA^iyj) = hA-iiyj) = ^{ijj) and h_{u) = h^^i^^) = (41) 

The quantity x + y + a is equivalent to ^ = —(1 + u!)~^k ■ {vi + V2) + a, and we 
find for the Isgur-Wise function 



1 



^ n2n n n n (^^) 



where Di = —2k -Vi — c. The normalization ^(1) = 1 is verified by setting vi = V2 = v 
and comparing with 
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Appendix B 

The first order correction terms to the weak decay form factors (defined in ( pH]) 
and (|21|)) may be expressed in terms of the universal functions as follows.]^ |13| 



7+e 



— + — ) (Xi + 2[1 - uj]x2 + 6x3) 



mc nibj \ 2 



A /uj-1 



2mc \uj + 1 
A /uj-1 



e + 2x1 - 4x; 



\u + 
A / 1 



A 



-[e-2e3] + 2xi+4[l-a;]x2 + 12x; 



[e + 6]+2x2 



2mc \uj + 
A 



2mt 



^—[{iu - l}e - 2^3] + 2x1 - 4x2 - 4x3) 
(e - 2^3 + 2x1 +4[l-a;]x2 + 12x3) 



In the notation of (|TT|) and Appendix A we obtain 

Sbhi 



a c (1 + uj) ■ {vi + f2) 
4 4^ DlD^D,D2 



X,; 



(43) 
(44) 



^(e + 2x1 - 4x3) + - 2^3 + 2X1 + 4[1 - u;]x2 + 12x3) (45) 



(46) 
(47) 



(4^ 



(49) 



' 4^^ DiD^D,D2 DlD^D.Dl 

a c Al + uy^k- {vi +V2) \ Yj 
4 4^ DlD^D,D2 J DlD^D,D2 



(50) 



where P or quantities are used in the expression for Sch^ according to whether the 
final state is pseudoscalar or vector, and Xi and Yi are given by 



i 








+ 


1 


1 l2 _ 
2 


ac 
4 


X+ 









ac 
4 


-x_ 


V 


1 










Ai 


1 


i+tj 


CO 

2(l+c^) 


y|^5 + Xa^ 


A2 










-2T + a(l + w)-iA;- (wi +^72) 


A3 


1 







-2t/ 



(51) 
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The quantities S, T and U are defined by 

On the other hand, proper normahzation of the vector meson "charge" graph at 
= gives in the notation of Appendix A 

n = 9R c 2ig-h)^ ^e-^ + cdy)^ a c k-v Yy 

^ A^DlD^D^^ DlD^D^ ^ ^ 2 ^ a' DlD^^D^ ^ DlD^D^ ^ ' 

where 

Yy = S-\k^-'^ (54) 

Setting vi = V = V2 va. (0) and (|5TD and comparing with (|53|) and (^4]) shows that 
^chAi = at u; = 1. 

An analogous argument using the normahzation of the pseudoscalar meson charge 
graph shows that (5fe/iAi = 0. Luke's theorem, XiiX) = X3(l) = 0, is therefore satisfied 
for arbitrary and with no extra constraints between model parameters. 

Returning to arbitrary uj we find that all the required relations which are supposed 
to hold between the 5bhi follow straightforwardly from the model. The following two 
relations which need to be checked for the 5chi are less trivial. 

Schy - 6chA, = TT-— — r (55) 

Z[L -\- id) 

ScihA, - hA, - hA, - 2(1 + uj)-^h^) = (56) 

(1 + ^) 

These are satisfied only if the following nontrivial identities are satisfied. 



2^ _ lj^2 _ k-(vx+V2) 

(57) 



c c 2(l+t>j) 

DlD^D,D2 



S+{l + uj){T + U)-p-{vi+V2) ^ c 

DIDMD2 2^ ^ 

We have confirmed these relations numerically. 

The functions X2 and ^3 are independent of g and h. We find 

_1 7 2 _ k-{vi+V2) 
C _ c 2(l+u)) 



DlD^D^D2 

C I C —?lL _l Q k-{v-i+V2) 

^ + ?3 , c ^ c l+o; 



(59) 



2(l + cu)+ DwId.dI • 
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